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1. Introduction 



The non-perturbative phenomena observed in the early 90's in the solvable models of 
non-critical string theories |l|-|7| and studied further in [|8|J| are now much better under- 
stood thanks to the recent regain of interest in this subject. In the papers [p!U|-p^|, the 



non-perturbative corrections to the string partition function were given a world sheet in- 
terpretation in terms of amplitudes of open strings attached to ZZ branes discovered by A. 
and Al. Zamolodchikov [|15| . The agreement between the matrix and CFT descriptions of 
the non-perturbative phenomena was observed not only the minimal (p, q) models of 2D 
quantum gravity but also the c = 1 string theory in the presence of vortex perturbations 
00 



which is believed to describe the 2D black hole [17 



As it was already pointed out in |L5| and then elaborated in []TT| , p!8| , |T9| , |20[| , the disc 
amplitudes on ZZ branes appeared to be the same as the disc amplitudes on FZZT branes 
taken at special complex values of the boundary cosmological coupling. A deeper geomet- 
rical understanding of this relation from the CFT side was achieved in the recent work of 
Seiberg and Shih []21| where the CFT description of the (p, q) models based on the ground 
ring structure [22,23]. It was shown that the ground ring relations lead to the same al- 
gebraic curve that appears in the matrix models approach. The algebraic curve found in 



21J1 gives the same representation of the correlation functions of the the FZZT brane in 



terms of Chebyshev polynomials as the one found in the loop gas models in [[M] and later 
in the two-matrix model in . The authors of gave a nice interpretation of the ZZ 
branes as degenerate (pinched) cycles of the complex curve. 

The disc amplitude on the ZZ brane associated with the degenerate cycle A mn is given 
by the contour integral of a certain holomorphic differential along the dual cycle B mn . The 
c = 1 version of this geometrical interpretation was proposed in |f26| . 

In the matrix approach the non-perturbative effects are produced by eigenvalue tun- 
neling amplitudes. Since the complex curve is determined by the shape of the effective 
potential, it is natural to seek a correspondence between the ZZ branes and the "non- 
minimal" saddle points associated with the local extrema of the effective potential. 

In this paper we show that the "non-minimal" saddle points are associated with the 
pinched cycles of the complex curve and express the instanton amplitudes through integrals 



along the dual cycles. Our results show that the geometrical picture found in [gTJ for the 
(p, q) critical points, actually holds for a general string background. 

We will consider the example of the two matrix model (2MM). It is known that all 
rational (p, q) models can be obtained by an appropriate tuning of the potentials in 2MM 



27,25]. The 2MM gives probably the most economical matrix description of all non-critical 



string theories with c < 1 matter content, such as pure gravity (c=0), exact solution of 
Ising on random dynamical graphs (c=l/2) [p8| , p9| or the 2D string theory on the self- 
dual radius [30-^2] • In case of a a polynomial potential, the planar limit of this model is 
described in terms of an algebraic (in general non-hyper-elliptic) curve [ |33| - |35 . 
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We evaluate the non-perturbative effects in the 2MM at the (p,p + 1) critical point 
by extending the quasi-classical analysis developed by F. David || for the 1MM. In order 
to obtain dimensionless quantities, we also calculate the normalized free energy for given 
complex curve using the Riemann bilinear identity. The results reproduce those obtained 
by Eynard and Zinn- Justin [Q from the string equation in the double scaling limit. 

We also consider the small perturbations of the algebraic curve around the (p,p + 1) 
critical point generated by order operators. We evaluate from the 2MM the leading order 
non-perturbative effects in the perturbed theory and compare them with the predictions 
of the boundary CFT which we extract from the known one point functions on the ZZ 
brane (normalized by the two point functions on the sphere to compare the dimensionless 
ratios). In the non-perturbed theory, there are four possible choices for the Liouville and 
matter boundary conditions, which lead to the same result. We found that this degeneracy 
is not lifted by perturbations by order operators, at least in linear order in the couplings. 
All four choices give to the same expression, which agrees with the one obtained in 2MM. 

Finally we discuss the c — > 1 limit and the instanton corrections in the Matrix Quan- 
tum Mechanics. We show that the instanton amplitudes in presence of vortex condensation, 
obtained previously using the equations of Toda hierarchy ]n,16| can be obtained as in- 
tegrals along compact cycles of the complex curve, very much as in the case of the c < 1 
theory. 



2. Spectral curve, free energy and instantons of the general 2MM 

In this section we will review the geometrical description of the planar limit of the 2MM 
|34|,|35| including the formulas for the free energy in terms of the integrals of a certain 
holomorphic differential along the closed A and B cycles on the algebraic curve of the 
model. We also will give expressions for the non-perturbative corrections in terms of 
similar integrals along non-trivial cycles passing through conical singularities, or double 
points, of the the algebraic curveR The conical singularities can be interpreted, following 
|21}1 , as degenerate, or "pinched", cycles of a curve of higher genus. These results will be 



used in the next section, where the (p, q) critical regimes of the 2MM will be considered. 

1 If the curve is defined by the equation F(x, y) = 0, then the double points are those for which 
dF(x,y) = 0. 
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2.1. Effective potential and saddle-point equations 
The partition function of the 2MM is defined as 



Z N = JdX dY e T>pcY-vCx)-v(Y)] = {1)n (21) 



where X, Y are the N x N hermitian matrices and 



q 



V(X) = J2 T k*\ V{Y) = J2 TkY k (2.2) 



k=l k=l 



are polynomial potentials. Using the representation of the 2MM in terms of the eigenvalues 



r N 

Zn = / II { dxk d W° eW-v^-vW) A N (x)A N (y) (2.3) 

k=l 



where Ajv(^) = Ylk>j=i( Zk ~ z i) ls ^ ne Vandermonde determinant, the ratio of partition 
functions Z^+i and Zn can be expressed as the double integral 



Zn+i 
Zn 



oo oo 

S c ft(x,y) 



dxdy e ~ b ^ x ' y \ (2.4) 



-oo — oo 



where S e ff(x, y) is the effective action for a pair of eigenvalues x, y: 

Ses(x, y) = -xy + V(x) + V(y) - log ( Det(x - X) Det(y — Y)) N . (2.5) 



The integrand in the r.h.s. of ( |2.4| ) is the expectation value of having one eigenvalue of 



the matrix X at position x and one eigenvalue of the matrix Y at position y. In the large 
N limit one can use the factorization properties ( Det A DetS ) = ( Det A ) ( Deti? ) and 
log ( Det A ) = ( Tr log A ), to write the effective action in the form 

S e «(x iy ) = -xy + $(x) + $(y), (2.6) 

with 

*(*) = V(x) - ( Tt log(x - X) ) N 
Hy) = V(y)-(Trlog(y-Y)) N 

and calculate the double integral ([2.4]) by the saddle point method. The saddle point 
equations 

x = $'(y), y = $'(x) (2.8) 
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determine not only the position of the saddle points, but also the potentials ( |2.7| ) them- 
selves. As follows from the studies of the two-matrix model by different techniques 
25| , ^4| , |37|j3^1 i, the solution of the model in the planar limit can be expressed in terms 



of a functional dependence between the complex variables x and y: 

x = X(y), y = Y(x) (2.9) 

where the functions X and Y are inverse to each other if considered as multivalued mero- 
morphic functions defined on their Riemann surfaces. On the physical sheets Y(x) — $'(x) 
and X(y) = $'(y) and therefore they satisfy (again on the physical sheets) the asymptotic 
relations 



Y(x) =V'(x) - N/x + o(l/x 2 ), x^oo 

_ (2.10) 

X(y) = V'(y)-N/y + o(l/y 2 ), y -> oo. 

These conditions and the fact that the two meromorphic functions ( |2.9| ) are inverse to each 
other and have no other poles except those at infinity determine them completely. 

2.2. Spectral curve 

The geometrical object behind is a complex curve F(x, y) = whose projections to 
the x and y complex planes are given, correspondingly, by the Riemann surfaces of the 
meromorphic functions y = Y(x) and x = X{y). For polynomial potentials ( |2.2| ), the curve 
is algebraic and its equation [|33|-|35| is a direct consequence of the asymptotics (|2.10[): 



F(x, y) = [y- V'(x)] [x - V\y)\ + P(x, y) = 0, (2.11) 
where P(x, y) is a polynomial of degree (q — 2,p — 2) 



V'(x] - V'(X) V'(v) - V'(Y) ■ 
P(*, y = Tr [X) V ^> V ^> UI1 - N . (2.12) 
x — A. y — Y 

The coefficients of the polynomial P(x, y) (the moduli of the spectral curve) are determined 
by the potentials ( |2.2| ) through the asymptotic conditions ( j2.10| ) as well as by the filling 



numbers Ni, N 2 , ... = N) associated with the cuts on the physical sheet [|34],j35 

As a real manifold the complex curve represents a two-dimensional surface of genus g < 
(p — l)(q — 1) — 1, with two punctures at x = oo and y = oo. 



2 A similar approach to the evaluation of the large ./V characters and heat kernels is used in 
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In the following we will consider the simplest situation when the algebraic curve has 
the topology of a sphere. We will also assume that 

q = p+l, (2.13) 

which is sufficient for describing the unitary non-critical string theories. Then the function 
( [2. 11|) has (p — l)(q — 1) — 1 double zeros, and the functions X(y) and Y(x) have a single 
cut on their physical sheets. 

The projection of the complex curve to the x complex plane, or the Riemann surface of 
Y(x), represents a branched p-covering of the Riemann sphere and thus splits the complex 
curve into p sheets. The physical sheet is the one that contains the puncture at x = oo. 
In addition, there are p — 2 cuts on the lower sheets that extends to infinity, where the 
Riemann surface has a branch-point of order p — 1. This branch point is the image of the 
puncture at y = oo, which we will denote by x = 6b. We will denote by Y^ k \ k = 1, 
the different branches of the function Y(x) keeping the label k = 1 for the physical sheet 
(Fig.l, right). 




Fig.l : Sheets of the Riemann surfaces of X{y) and Y(x) for the one-cut solution 

with p = 3, q = 4. 



In a similar way the Riemann surface of X(y) splits the punctured sphere into q sheets 
X^\ ...,XW (Fig.l, left). Then the spectral curve can be rewritten as |34j] 



F(x, y)~f[(y- Y (k \x)) ~ f[ (x - XW{yj) . (2.14) 



k=l i=l 
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We will refer to the points x = oo and x = 6b as the north and south poles of the 
sphere. The punctured sphere is characterized by the two cycles A and B dual to each 
other. The cycle A goes along the equator and the cycle B connects the north and the 
south poles. We assume that the sheets can be labeled so that there is exactly one cut 
connecting the /c-th and the k + 1-th sheet. The first and the last sheet of the Riemann 
surfaces of Y(x) and X(y) have one cut and the other sheets have two cuts. 

There is always possible to find an uniformization parameter 10 that globally 
parametrizes the complex curve [Z7],[25fl, 



k=~i j=~i 

One of the coefficients can be chosen arbitrarily because of the symmetry with respect to 
rescalings of u. The rest p + q + 1 coefficients are determined through ( |2.1U| ) as functions 
of the p + q couplings Ti, T p , T\, T q in ( |2.2| ) and the number N of eigenvalues. The 
punctured sphere is parametrized by the complex plane u. The north (south) pole of the 
sphere then corresponds to the point u = oo (u = 0). 

2.3. Perturbative free energy in terms of the spectral curve 
The saddle-point equations ( |2.8| ) mean that for some k and / 

x = X {k) (y), y = Y^(x). (2.16) 

Inverting the second equation we find that for some j ^ k: X^ k \y) = X^\y) which is 
satisfied when y is at one of the endpoints of the physical cut of X(y). The same is true 
for the variable x. The "perturbative" saddle points are at the endpoints of the physical 
cuts of the functions X(y) and Y(x) i. 

By the asymptotics ( |2.10| ), the integral along the A-cycle is equal to the number of 



the eigenvalues: 

N = (b ydx. (2.17) 



A 



As was shown in |J(and generalized to the filling not only of the maxima but of of all 
extrema of the potential ||44J| ) for the one-matrix model and then for the two-matrix model 



The same algebraic curve can describe several matrix models; they correspond to different 



real sections that define different sets of local minima of the effective potential [35,32]. Note that 
in the Normal matrix model, which has the same complex curve as the Hermitian two-matrix 
model, the saddle point is described by a closed contour and not by isolated points. This is 
possible because the Normal matrix model is described by different real section of the complex 
curve [ 43 1 . 
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in [35|], the integral of the function Y(x) along the S-cycle is equal to the derivative in N 
of the planar contribution JF to the free energy 



F = \ogZ N . (2.18) 

Indeed, the leading contribution to the integral (^]) is given by the saddle-point value 
of the effective potential ( |2.6| ). Let x',y' be related by y' = Y^(x') and x' = X^(y'). 
There is always such a point along the cycle B. Then we write 

d N ^o = -S e s(x',y') 




(2.19) 



The formulas (|2.17[ ) and ( |2.19| ) determine the free energy in the planar limit in terms 
of the two main cycles of the punctured sphere. These formulas are geometrical in the 
sense that they do not depend on the choice of the coordinate patches on the complex 
curve. In the next subsection we will obtain similar formulas for the non-perturbative 
instanton contributions. 



2.4- Leading order non-perturbative corrections 

Besides the perturbative saddle point there are also other saddle point solutions, which 
describe the non-perturbative corrections to the free energy [|5],|6|||. The meaning of the 
non-perturbative corrections depends on the physical context. Typically they describe the 
decay of a metastable ground state caused by tunnelings of eigenvalues under a maximum 
of the effective potential (eigenvalue instantons) . 

A comprehensive description of the instanton effects in the one-matrix model has 
been done in ||. Once the effective action is known, one can generalize the analysis of 
H to the case of the two-matrix model. Here we will restrict ourselves to the leading 
non-perturbative corrections, which allow geometrical description in terms of the spectral 
curve. 

The "non-perturbative" saddle points are double points of the complex curve repre- 
sented by the pairs x = Xku V = Vki such that 

y = Y {k \x), x = X {l) {y). (2.20) 
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At these points dF(x, y) = and the effective action ( pT6|) has vanishing derivatives in x and 
y. The pairs (1, 1) and (fe, 1) should be excluded because they describe the same point of the 
curve (the second function is the inverse of the first). The pair (2, 2) should be excluded as 
well because it determines perturbative saddle point. The number of the remaining pairs 
(k, I) is equal to the maximal genus of the complex curve, <7 max = (p — 1) (q — 1) — 1. Among 
these saddle points there are (p — l)(q — l)/2 maxima and (p — l)(q — l)/2 — 1 minima 
of the effective action. In the last case the integration contour should be distorted in the 
complex plane as explained, say, in Section 2 of p5| . 

Each pair (xki,ykl) correspond to two different points of the complex curve which 
can also be considered as a remnant of a collapsed handle, or vanishing cycle Am of the 
complex curve in generic position |2l]. Such a vanishing cycle is associated with a pair of 
coinciding branch points of the Riemann surface of Y(x) or X(y). Let us denote by Bki 
the corresponding S-cycle, which connects the two poles of the sphere through the pinched 
cycle Aki- Proceeding as in (|2.19| ) we can write the effective action at the (k,l)-th saddle 
point 

S e fi = /_ ydx = -d N T Q + S kh (2.21) 



where 

S k i = I ydx. (2.22) 



is the effective action associated with the compact cycle Bki passing through the pinched 
cycle and the cycle A (Fig. 2). 




Fig. 2 : The only cycle B23 for pure gravity (p = 2, q = 3). It can be viewed either 
as an open contour connecting the two double points on the sphere (left) or 
as a closed contour going through a pinched cycle of a torus (right). 
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If we take into account both minimal and non-minimal saddle points, eqs. (2T3) and ( |2.21| ) 
give 

2, 



Z 



N 



c ° \1 + J2 c kie- Skl (2.23) 



k.i 



where the sum is taken over all pairs (k, /) discussed above. Written for the free energy 
T = logZ/v, this formula reads 

d N F = d N F + c kl e~ Sk '+... (2.24) 

(fe,/)#(2,2) 

As was explained in ||,0, the non-perturbative corrections should be understood not as 
an improvement of the 1/N expansion but rather as a difference between the two free 
energies corresponding to the two ways of distortion of the integration contours in the 
matrix integral into the complex plane. 

We believe that in any matrix model described by algebraic curve the non-perturbative 
corrections are given in the planar approximation by the formula fl2.24f )il. In the next 
section we will apply (|2.24| ) to the critical regimes of the 2MM, in which case it reduces to 



the similar formula obtained in 121 



3. Explicit results for the (p, p + 1) critical points 

3.1. The scaling limit near the (p, q = p + 1) critical point 

The noncritical string theories are described by the critical points of the complex curve. A 
critical point can occur when a branch point comes close to a double point or to another 
branch point. Here we will examine the "maximal" critical point that arises when the right 
branch point point of the physical cut of Y(x) coalesce with the p — 2 branch points on 
the lower sheets. At this point the equation of the curve ( |2.15| ) takes the form [^5j 

X(u) = N c ( 1 ~ a; ) P > Y{u)=N c u(l--) . (3.1) 

In the vicinity of the origin the complex curve degenerates to 

y p = x q . (3.2) 



The solution ( |3.1| ) corresponds to certain choice of the coupling constants in ( |2.2| ), all 
of order of iV c . One can introduce a one-parameter deformation of this singularity by 
changing the number of eigenvalues to N < N c . The difference N c — N is proportional 



We have checked that ( 2.24 ) holds also for the ADE matrix chains defined in [46|. 
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to the cosmological constant \i of the corresponding (p, q) string theory. By introducing a 
small cutoff parameter a and rescale the variables as 



x y N - N c 

► > y, — : — — 

ai aP aP+iN r 



9s 



(3.3) 



we blow up the the vicinity of the (deformed) critical point so that all the cuts become 
semi- infinite (Fig. 3). 





/ / ** '*'*, \ 
/ * * \ 
f i OO \ ] 










..... x il} \ 




Fig. 3 : Sheets near the (p, q) critical point with p = 3, q = 4, in the 
parametrization (|3.1|). On the right, the blown up scaling domain. 



The physical cuts of the functions y = Y(x) and x = X(y) extend to infinity along the 
intervals — oo < x < —2M and — oo < y < — 2M respectively, where M ~ /U 1 / 2 and 
M ~ M q l 2v . We will normalize X and Y so that M, M and \i are related by 



(3.4) 



A very useful parametrization in the scaling limit is given by expanding X and Y in the 
Chebyshev polynomials of a third variable 

z = 2£ cosh 9, 

which appears naturally in the formalism using the dispersionless KP hierarchy ||25|| : 

x = 2T p (z/2) = 2£ p coshp#, 
y = 2T q (z/2) = 2£ 9 coshg#. 



(3.5) 
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To avoid the subtleties that arise in non-unitary theories, we will assume below that 
p = q + 1. The parametrization (|3.5| ) unfolds the branch points of the functions Y(x) 
and X(y) and allows us to work with entire functions of 9. The sheet structure is shown 
in Fig. 4. The /c-th sheet of the function Y(x) is parametrized by the semi-infinite strip 

k k — 1 

Re0>O, -7r<|Im0|< 7r 

p p 

and the Z-th sheet of the function X(y) is parametrized by the semi-infinite strip 

Re#>0, — - — 7r <C I Im (tv — $)| <C — — — — 7r. 

p+1 ' V n p+1 

At the critical point the parameter z is related to the parameter u in by 

iv = 1 + az. 



Y (2) 



271 v' jJ 



ImG 



X 



(4) 



X 



(3) 



X 



(2) 



X 



Re9 



Y <2> 
Y (1) 



X 



(2) 



X 



(3) 



X 



(4) 



Fig. 4 : Sheets in the 6 parametrization for p = 3, q = 4. 
All functions are symmetric under reflection 9^—6. 

In any (p,p+ 1) theory of 2D gravity with matter central charge c = 1 — , 6 +1 y , the function 
y(a;) defined by ( [3.5|) gives the loop amplitude, that is the disc amplitude with a marked point 
on the boundary, with the boundary cosmological constant x and bulk cosmological constant /U. 
The function X{y) does the same for the dual (p + theoryll. The expression (|3.5| ) has been 



The duality among (p, q) and (q,p) theories was first observed in the matrix models in [47,48| 
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first derived in p4}| in the context of the ADE loop gas models. Its interpretation from the point 



of view of Liouville theory is given in [49| 



The two potentials ( |2.7j ) can be calculated in the scaling limit by integrating ( p,5|) . This 
gives (q = p + 1) 



$ = 2p £ 



p+q 



$ = 2g £ 



p+q 



cosh(p + q)9 cosh(q — p)9 

p+q Q—P 
cosh(p + q)9 cosh(g — p)0 



(3.6) 



p+q q— p 

and one can check that the effective action is constant on the complex curve (equal to zero in our 
conventions): 

$(s) + 4>(y) = xy. (3.7) 
To make connection with the boundary Liouville theory we introduce the parameters 

b = \/ p/q, itbs = ^fpq 6. (3.8) 

Then the equation of the complex curve becomes 

x = cosh rrbs, y = 2/x 1 ^ 2b cosh-7rs/6. (3.9) 

3.2. Perturbations by relevant operators around the (p,q = p + 1) critical point 

A generic perturbation of the (p, p+ 1) critical point is described by a curve that behaves as ( fj.2| ) 
at infinity. Such a curve has the parametric form 

x = x(z) = z p + Xp-iz 9 ^ 1 + x p -2Z p ~ 2 + . . . + y , 

(3.10) 

( \ P+ 1 i P— 1 i P-2 ,i V ' 

y = y{z)=z y +y p -!Z y +y P -2Z F +... + y . 

The curve ( |3.10| ) has only one singular point at infinity where the function y(x) has a branch 
point of order p and the function x(y) has a branch point of order q = p + 1. 

The original Toda integrable structure of the two-matrix model is characterized by two 
singular points, oo and 6b, and the spectral curve is determined by two asymptotic conditions 
( 2.10D associated with the two punctures. For degenerate curves of the form ( [3. 10] ), the relevant 



integrable structure is that of the p- reduced dispersionless KP hierarchy pOJ , which has only one 
singular point y = x = oo. Therefore in the scaling limit it is sufficient to introduce a single local 
coordinate in the neighborhood of the puncture, say x, and the corresponding set of couplings 
t = {tn}^! defined by the coefficients of the powers x n ^ p ~ 1 , n = 1, 2, in the Laurent series of 
y{x) at infinity: 

y(x) = d x ^ = -J2( n ^ x^ 1 " 1 + vn x^- 1 ) . (3.11) 
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The curve ( 3.10 ) is described in terms of the non-zero couplings t\, ...,t p -i,t p+ \, ...,t2 P +i- The 
coefficients v n in ( [3.11 ) are functions of these couplings. The integer powers in x do not have a 
branch point at infinity and therefore the sum is restricted ton/0 (mod p). 

The integrable perturbations associated with the couplings t n are generated by the Hamil- 



tonian flows of the p- reduced KP hierarchy. Following Krichever [ 50 ] , one can associate with the 
Hamiltonian coupled to t n a meromorphic differential dH n (x,t). The form of the classical Hamil- 
tonians are determined by the complex curve ( |3.10 ). The differential of the effective potential 
&(x, t) is given by 



d<£ = ydx + H m (x)dt m . 



(3.12) 



The dispersionless KP hierarchy can be interpreted as a classical Hamiltonian system with Poisson 
bracket 



{/,<?} 



df dg df dg 



(3.13) 



dz dt\ dt\ dz 

The differential (gl2|) then can be considered as a generating function for the canonical trans- 
formation between the phase-space coordinates (z,ti) and (x,y). Given the curve (13.10|) and the 
expansion ( p.ll| ) at infinity, the expressions for the classical Hamiltonians are 



H„ 



[X 1 



i /p] 



(3.14) 



where [ ]+ denotes the non- negative part of the Laurent series of the function X(z) at infinity. 
It follows from (3.12) that the parameter z in the definition of the curve ( |3.10 ) can be identified 
with the first Hamiltonian Hi: 



Hi = [X 1/p ] 



(3.15) 



It is technically convenient to use the parametrization ( |3.10|) and consider z = Hi as a global 
coordinate on the complex curve. From ( |3.12| ) and ( |3.14| ) one finds the expression of the effective 
potential 

$ = ^ ] H m (z,t) t m - ^3 ^Q-j 



i>i 



Given the relevant couplings ti, ...,t p -i,t p +i, ...,t2 P -i, from ( pM6| ) and ( 3.14 ) one can then eval- 
uate the 2p — 2 coefficients defining the curve ( |3.10 ). 

Now let us return to the special case of the curve ( |3.5| ). The only deformation parameter 
here, the cosmological constant /j,, can be identified, up to a normalization, with the coupling t\. 
Plugging ( [3.5[ ) into (3.14) one finds for the classical Hamiltonians [25] 



H m = 2£ m cosh m9, m = 1, 2, 2p — 1; 



H 2p+1 = 2f p+l ^cosh(2p + 1)9 + cosh 6^J 



(3.17) 
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Comparing with the expression ( |3.16[ ), we conclude that the curve (|3.5| ) describes the point 

P 

t2p+l = , tx = —(p+l)fjL, Others = (3.18) 

ZJ) -f- 1 

in the space of couplings. The relation between /x = £ 2p and t\ follows also directly from the 
classical string equation {x, y} = 1 applied to the solution fl3.5D - 



3. 3. Free energy and two-point functions 



Here we will calculate the two-point functions of relevant operators on the sphere for the back- 
ground ( |3.5| ). The pair of dual variables ji and d^T are expressed in terms of the complex curve 
through the formulas ( |2,19| ) and ( |2.17| ), In terms of the rescaled variables defined by ( |3.3| ) these 
formulas read 



and 



-d a T 



ydx 



d$ 



B 



ydx 



d$. 



(3.19) 



(3.20) 



The general formulas for the 2MM with polynomial potentials, of the type obtained in [51] 
cannot be immediately applied for this task in the critical limit. In the scaling limit the integrals 



in ( 3.19 ) and ( 3.20 ) diverge and need regularization. One could do it explicitly by making the 
cut finite. It is however possible to extract the necessary information from ( |3. 19 ) and ( |3.20 ) 
without any regularization by using the Riemann bilinear identities (RBI). RBI follow from the 
fact that there are no meromorphic (2, 0)-forms on one-dimensional complex curves. For each pair 
of holomorphic differentials dfli and d^2, 



= / rffiiA dft 2 



(b d£li (b dQ,2 — (D dQ,2 (p dfli — res (^1^2) - 

J A J B J A J B 



(3.21) 



In particular, when 

dQ-i = d tl Y(x) dx = ddfj,Q, dQ 2 = d m d n Y(x) dx = dd m d n Q, 
equations ( |3. 19 ) and ( 3.20| ) yield 



(3.22) 



— (p dd^(x) ■ d m d n $>{x) = d^dmdnF. 



(3.23) 



A rigorous derivation of this identity can be done using the the fact that the partition function 
of the two-matrix model is a r-functions for the KP integrable hierarchy, which means that the 
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coefficients v n in ( |3.11| ) are related to the free energy by v n = d n T . This implies the following 
identities, which have been derived in [p0| (for more explicit derivation see [52],[ 



dmdnJ 7 : 



2ni 



2tu 



x n/p d<f>. 



x n/p dH n 



(3.24) 



dkdmdnT = -7^— f x n/p dd k H m = — !- <b (d k H m ) x dH n . 



Since \l ~ ti, equation ( 3,23| ) is a particular case of the last of these relations. 

Our aim is to calculate the two-point function d m d n J- in the classical background ( |3.5| 
where the classical Hamiltonians are given by (3.17). Inserting 



(a u \ ^-2 P +m smh(p - m)9 n . 



into ( |3.24| ) we get 



d^d m d n T = 2mn £ T 



n -2 P f d6 s'mhnO sinh(p — m)9 



27rz psinh p9 

We evaluate the integral along the contour lm8 6 [0,2-7r), He 9 — > 00, which gives 

do t) 77 TtlTl J2{yi — p) ^ 

P 

Integrating with respect to \x = £ 2p we find for the two-point function 



(3.25) 



(3.26) 



(3.27) 



d 2 n T = n£ 



2n 



(3.28) 



For k = 1 we obtain from here for the string susceptibility 



(3.29) 



3.4- Eigenvalue instantons at the {p,p + 1) critical point 

As was discussed in section 2.3, the eigenvalue instantons are associated with the the points ( |2.20 ) 
of the curve where conical singularities occur. These points are given by the non-trivial solutions 
(9^9') of the equations 

x(9) = x(9'), y(9)=y(9'), (3.30) 
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where the functions x{9) and y{9) are defined by ( |3.5[ ). The solutions of (3.30) are 9 = 9 mn , 9' 
9 m -n (m = 1, ...,p - 1; n = 1, with 



ITT 



rn 
V 



n 



p + 1 



(3.31) 



The equations (|3.30|) can be written also in the form ( 2.20 ), with k = [\m — n^j\] and I = 
[{m^^- + n\]. Due to the symmetries 9 — > —9 and 9 — > 9 + 27r, the points with parameters 9 mn 
and 0p-m,p+i-n describe the same double point of the complex curve. Therefore in the scaling 



(p-l)(g-l) 

action. The (p ~ 1) n (g - 1) 



limit there are only KV L) ^ q — — double points which correspond to local maxima of the effective 



1 double points that correspond to local minima are sent to infinity after 
the blow-up. In Fig. 5 we give the plots of the closed contours in the (x, y)-plane for p = 2, 3, 4. 
Some of the cycles have backtracking parts, which can be deleted. The effective action associated 
with each cycle is proportional to its (algebraic) area. 




Pure gravity (p = 2) 




Ising model coupled to gravity (p = 3) 




Tricritical Ising model coupled to gravity (p = 4) 
Fig. 5: The instanton cycles for the (p,p+ 1) non-critical strings with p=2,3,4. 



The instanton contribution to the effective action corresponds to the non-trivial cycle B mn and is 
given by ( ggg ) 

S mn = I Y(x)dx = $(e rnn )-<f>(0 m ,-n). (3.32) 



Evaluating the difference with $ given by (|3.6|), we get 



S mn = £ 2p+1 sin(7rm/p) sin(7rn/g). 



(3.33) 
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We can also calculate the first order corrections^ 5S mn = ^2, k tkdkS mn to the instanton 
effective action in presence of perturbations by the order operators Ok, k = 1, ...,p — 1: 

dkS mn = H k {9 m ,n) - H k {9 m - n ) = -4£ fe sin sin "JLJU. (3.34) 

P 1 

To compare with the world sheet CFT we need to prepare a dimensionless ratio that is not 
sensible to the normalization of x, y and fi. We introduce the dimensionless ratio of ( |3.34 ) and 

(EH) 

(it) dkSmn 4 . -Kkn . irkm . orN 

r m ' = — = F sin sin 3.35 

Vk P q K ' 

For k = 1 the result coincides with those found in this particular case in H and reproduced within 
the Liouville CFT approach in [O] . We will show now that this result can be reproduced also for 
any k from the CFT. 



4. One-point functions on ZZ branes and instantons 



Let us compare the resu lt ( 3.35| ) to the normalized one-point functions on ZZ branes in a Liouville 
theory with Q = yj £±I -|- ^/_£_ ) i n the same spirit as it was done in ]TT| , |T4| ] . The Hamiltonians 
Hk are described in the CFT approach by the product of matter and Liouville operators 



Hk 



O k ■ e 2ak * 



(4.1) 



Dis 



We assume the most general boundary conditions for the matter and Liouville fields. The matter 
boundary conditions are labeled by the entries of the Kac table, which we denote by (m, n) ~ 



(p — m,p + 1 — n) with 1 < m < p — 1 and 1 < m < p [54]. The ZZ boundary conditions for the 
Liouville field (m',n') with m',n' > 1 EI]]. We will need the following formulas: 

- The structure constant for the disc one point function of the matter order operator Ok in 
presence of the boundary condition m,n [55|,56|: 



(Ok) 



matter 
ra.n 



pip + 1) 



1/4 gin Trmfe gin Trnfe 
V P+l 



sin sin ^- 

p p+i 



(4.2) 



The disc partition function in presence of m',n' ZZ-boundary condition of the Liouville vertex 



operator e 2afe *, a.k 



t • oin Kkm' ■ Trkn' 

Liouv sin sin ■ 



P+l 



sin — sin 

p p+i 



6 The pre-factors Cki in (2.23) were conjectured in fi"4[| to behave as c r , 
theories, as is known from some explicit examples (see 0,^4|). 



(4.3) 

~ gV 2 for all c < 1 
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where 



LlOt 



- M k/p - 

7T 



2 3 / 4 7Tfc 



V^+i)r(i + i A_)r(i + |) 

and the constant M is related to the Liouville bulk cosmological constant \xl by 



(4.4) 



M 



p 



p + 1 



7(35) 



r(x) 
rri-x) 



(4.5) 



-The Liouville two-point function on the sphere [57,5q| (see the Appendix and the eq. (3.45) of 
14] for this particular formula): 



sphere 



*Vp(p + i) 

2ttp 2 



M 



2fe/p 



7 1 



P+1 



7 



(4.6) 



Now we have to see which matter and Liouville boundary conditions can reproduce the 



dimensionless ratio ( 3.35| ) obtained from the two-matrix model. The explicit formula for 



obtained by combining (|4.2j ), ( |4. 3|) and ( |4.4j) is 



m,n:m ,m 



matter , a A \ 

\® k 'm,n \ e )m'n' . TTmk . 7mk . TTm' k . TVIl' k 

' — = pk sin sin — —r sin sin 



with 



pk 



(e^e^) sphere 

Matter , , , Liouv 

(Mia •< eQfe ') 11 



p + i 



v 



p + i 



l e a k ct> e a k 4>\ 101 
\ / sph> 



Liouv 
ere 



a/2A- 



r (1 + A) r (1 + fc/ri f-^4i^7 (1 - pfr) 7 (-Vp) 



1/2 



(4.7) 



(4- 



4 

We see that the matrix model result matches with a particular subset of boundary conditions. 
The exponent of the tunneling amplitude in the matrix model associated with the cycle m, n is 
reproduced by r„ , , with any of of the following four choices 



m,n;m,n 



m, n; 1, 1 
1, 1; m, n 
1, n; m, 1 
m, 1, 1, n 



(1) 
(2) 
(3) 
(4) 



(4.9) 



7 See the discussion on the origin of the first factor at the end of subsection 2.2 of the paper 
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In the case of a non-perturbed theory (tk = 0), the first choice was considered in []l2|,[l4j and the 
last two were considered by Martinec 

Our first-order calculation shows that this degeneracy is not lifted by perturbations by order 
operators. The matrix model does not distinguish between the four choices in ( |4.9| ) and there is 
no convincing physical argument that allows to single out one of them. 

It has been argued by the authors of [21 that in the (p, q) string theory the \m, n; 1, 1) branes 



with 1 < m < p — 1, 1 < n < q — 1 and mq — np>0 form a complete set of distinct physical states 
with ZZ-type boundary conditions, which they called "principal branes". The other branes should 
be thought of as multi-brane states formed out of these elementary ZZ branes. The degeneracy 
we observed suggests that all four choices in ( f4.9|) provide a possible basis of "principal" branes 
and each of the choices describes the same set of physical states. 

This statement seems less strange if we recall the following two facts. First, it follows from 
the exact expressions for wave functions of the FZZT and ZZ boundary states [5£,15|, that the 
ZZ state (m, n| can be obtained as a difference of two FZZT states {9\ 

(m,n\ ~ (9 m ,n\ - (0rn,-n\ = (#m,n| ~ (0- m ,n\, ( 4 -10) 

where the angle 9 is related to the boundary cosmological constant as ^lb ~ y^cosh^^) and 9m,n 
is defined by ( |3 .31 ). This fact has been further explored in [6^,19,31,11]. On the other hand, it 



is also known that the target space dimension in the rational string theories is associated with 
the imaginary direction of the uniformization parameter 9. Namely the points of the discrete 
target space are labeled by the p — 1 cuts of the Riemann surface of the function Y(x) or the 
q — 1 cuts of the Riemann surface of the function X(y). The discontinuity of Y(x) along the 
m-th cut is given by -2irip {m) (9) = Y(9 + m^) - Y{9 - inf). From the point of view of CFT, 
the m-th cut describes a matter boundary condition of type (m, 1). Similarly, the discontinuity 
of X(y) along the n-th cut is given by -2-Kip^ n \9) = X(9 + iiv^) - X{9 - m~) and describes 
the matter boundary condition (l,n). The matter boundary conditions (m,n) has never been 
studied from the matrix point of view, but it is plausible that they can be described in terms 
of the disc partition function &(9) with boundary parameter shifted in the imaginary direction 
9^9 + 9± m ,± n - Thus it seems that the translations of the boundary parameter by 9± m ,± n can 
be interpreted either as projection to (m, n) TIL boundary condition for the Liouville field, or as 
projection to (m, n) boundary condition for the matter field. Of course this statement needs to 
be better understood. 
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5. The c — > 1 limit and comparison with Matrix Quantum Mechanics 

The c = 1 limit is obtained by taking b 2 = -» 1. In this limit it is convenient to use the 
parameters s and M related to 9 and £ as 

7rbs=p6, M = i p . (5.1) 

The complex curve ( |3.5| ) 

x(s) = Mcoshnbs, y(s) = M l ^ b coshns/b (5-2) 

becomes in the limit p — » oo 

x(s) = Mcosh7rs, y(s) = x + —\x log M + ttsM smb. -its]. (5-3) 

P 

In this limit it is convenient to redefine the variable y by subtracting the linear in x term and 
rescaling by 1/p: 

y(s) = irsM sinh ns. (5-4) 



The function (5^4) gives the continuum limit of the resolvent in the ADE matrix models [62], 



which describe particular sectors of the c = 1 string theory. 



In the description based on the Matrix Quantum Mechanics |63|, the natural variable is the 



canonical momentum p conjugated to the eigenvalue x, which is related to the resolvent by 



*«) = »('"') ^(°-"'> (5.5) 



or, in terms of the parameter s, 



p ^ = y(s + i) y(s — = Mginh ( 7rs ) ) x = Mcosh(Trs). (5.6) 

The relation 

p 2 — x 2 = /x, /j, = M 2 , (5.7) 

is the equation for the classical phase space trajectory in MQM. In the world-sheet CFT this 
relation appears in the context of the Witten's ground ring 

In the c = 1 theory it is more natural to consider the chiral variables 

x+=p + x, x-=p — x (5-8) 

which describe the left and right moving tachyons. The variables x± can be considered as a pair 
of local coordinates covering the two-dimensional sphere. The punctures at the north and the 
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south poles of the sphere correspond to the points x+ = oo and X- = oo and the two charts are 
related by the equation of the complex curve Q5.7| ) 



X+X- 



(5.9) 



As in the case of the (p, q) string theory, the perturbations are introduced as the asymptotics at 
the two punctures. The complex curve deformed by tachyon sources ^2 n>Q t± n (x±) R is found in 
[IJ. The analog of is 

d<E> + = x-dxA 



x+dx- 



^ H n dtn 
n>0 

^ Hfi dtn , 

n<0 



(5.10) 



where the Hamiltonian H± n is given by the positive/negative part of the Laurent expansion of 
the function x±{uj), uj = e 7 ", plus half the constant term. Note that here we should keep both 
expansions since the singular points x+ = oo and X- = oo are distinct. The parameter oj is 
expressed as a function of x+ or X- as u) = e l9f '*+^+- ) 
(l3T20| ) turn to 



e ia M*-(*-). The conditions (|3T9|) and 



x-dxj 



B 



x-dx+ 



(5.11) 



d$~ 



where the compact cycle A goes around the equator and the non-compact cycle B connects the 



two punctures [32]. 

Now let us concentrate on the instanton corrections in MQM compactified at radius R and 
their CFT interpretation. In the non-perturbed theory described by (]5.9| ) the non-perturbative 
corrections follow from the integral representation of the free energy 



n = l 77 = 1 



-2 TV n. R \i, 



(5.12) 



The world-sheet description of the two kinds of exponential terms was discussed in |14[ | . The terms 
e -27rn/x correS p 0n d t Dirichlet boundary conditions for the matter field. Due to the translational 
symmetry of the target space, the instantons are labeled by only one number n (the other one 
is redundant). The integral ( |2.22 ) is replaced by an integral along a cycle going n times around 
the neck of the hyperboloid. Thus the n-instanton solution is associated with a contour winding 
n times around the neck. The original algebraic curve is actually the universal covering of the 
hyperboloid, due to the logarithmic dependence of y(x). The curve wraps infinitely many times 
the hyperboloid x + x~ = [i. Similarly, the terms e ^ 2nn ^ R correspond to Neumann boundary 
conditions for the matter field and correspond to the cycles of the dual curve describing the 
vortex excitations. It was argued in [16] that the Liouville boundary conditions for the n-th 
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instanton correction are given by the (l,n) ZZ boundary state. The argument of fl6|] was based 
on the comparison of the first order perturbations in MQM and the world sheet CFT. 



Here we will show that the instanton corrections found in [14,16] are related to the disc 



partition function by a formula similar to ( 3.32 ). This statement is trivial when the theory is not 
perturbed. Indeed, in this case 



and 



1/2 ±1 1/2 ±7TS 

X± = fj, LO = U e , 

= ^logO~ 1/2 a;+) = fins - -ulog 



$+(s = in) — <£+(s = —inn) = 2irinu, 



(5.13) 



(5.14) 



which reproduces the first kind of instanton corrections. By construction the potential <&±(x±) are 
the generating functions of the tachyon operators and as such describe Dirichlet branes with fixed 
time position. The second kind of instanton corrections is reproduced by the effective potential 
3>± in the dual theory, which describes vortex excitations and therefore Neumann branes wrapping 
the time circle. 

Below we will consider the effect of the tachyon perturbations with t% = t-i = A, tothers = 0, 



which affect only the Dirichlet branes. Then the instanton corrections take the form [14,16] 



— 2irnR^t 



with 



1 A 2 1 

/n(/i,A) = 27m + 4Asin(7m/i?) / uA- 1 + sin(27rn/i?) //^~ 2 + ... (5.15) 



The exponent ( |5.15| ) can be presented again as 



n fn{n, X) 



<£> + (s = in) - $+(s = -in) 



(5.16) 



where 



$+(>+) = ^log(^ 



-1/2, 



2A/x« (x + ) R 



A 2 A 



2R 2 



(x + )-k + ... 



(5.17) 



is the effective potential in presence of perturbation. It can be formally interpreted as the disc 
partition function with a "crural" (and therefore non-local) boundary condition labeled by the 
value of x+. The relation ( 5.1 6| ) can be generalized to the case of finite perturbation [26]. 
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6. Discussion 

In this work we derived the non-perturbative corrections in the unitary (p,p + 1) models 
of 2D quantum gravity from the collective field theory for the two-matrix model. We used the 
approach developed in |5]|| for the one-matrix model. The results have a nice algebro-geometric 
interpretation: the contribution of an instanton corresponds to the integrals over the cycles start- 
ing and ending at two different sheets of the Riemann surface of the algebraic curve of the model, 
at the points corresponding to a pinched cycle connecting two sheets. In the CFT description this 



interpretation was proposed in [21| where the curve appeared from the ground ring relations. 



We also compared the new results with the Liouville CFT (generalizing the successful com- 



parison done in |14| (see also [11] where the idea of such a comparison was proposed). The CFT 
interpretation of our matrix instanton calculation is the following: we computed the one point 
functions of primary fields on the ZZ brane already found on the CFT side in [15]. As usual, 



the normalizations of operators are different in matrix and CFT approaches, and it only makes 
sense to compare the dimensionless ratios of various quantities. For that we had to calculate the 
free energy and the two point functions of the primary fields on the sphere in the same instanton 
approach. This calculation appears to be subtle due to the singular nature of the critical algebraic 
curve. We succeeded to do it using the geometry relations [p|,p,44] defining the free energy of the 



model through holomorphic integrals over the curve, combined with the bilinear Riemann identity. 

Our method may be not as straightforward as the old method based on the string (KP) 
equations but it is more transparent geometrically and might have a wider range of applications, 
especially for the matrix models where the method of orthogonal polynomials does not exist but 
the curve is known. 

The instanton method also clarifies the geometrical meaning of the relation between FZZT 



brane and ZZ brains, as mentioned in [pi,19|: they correspond to different choices of the contour 
in the same holomorphic integral: for the FZZT brane the contour starts at any point on the curve 
(corresponding to the complex boundary cosmological constant) and goes to infinity, whether as 
for the ZZ brane the contour connects two different critical points corresponding to the pinched 
cycle of the algebraic curve. 

Our results show that the agreement between the matrix model and CFT calculations es- 



tablished in [14] for the (p, p + 1) critical points, holds also in presence of perturbations by order 
operators, at least in the linear order in the couplings. The most interesting physical outcome of 
our calculation is that the results do not depend on the choice of the ZZ branes: (m, n), (m, 1), 
(1, n) and (1,1) ZZ branes give the same instanton effects being combined with (1,1), (l,n), 
(m, 1) and (m, n) matter branes, respectively. This degeneracy is needed for the self-consistency 
of the matrix interpretation of the ZZ branes as non-perturbative effects since the difference in the 
results would mean that we were missing the matrix model description of some of these branes. 
This degeneracy is still to be understood within the Liouville CFT of 2D gravity. 

We also applied our method to the study of the instanton effects the c = 1 string compactified 
on an arbitrary radius and perturbed by relevant vortex operator. We showed that the instantons 
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have here essentially the same geometrical meaning in terms of the holomorphic integrals along 
the algebraic curve as for c < 1 models. We also identified the multi-instanton contributions in 
this case as multiple windings of the integration contour along the cycles of the curve. It is known 



that this model gives the black hole realization proposed in [17]. The non-perturbative corrections 
play there the crucial role for the understanding of black hole physics in the near horizon strongly 
coupled area. Our method opens the way for the geometrical instanton interpretation of the 



energy of the ZZ-brane in the black hole calculated in [14] by the matrix approach using Toda 
equations of [17]. 

An interesting continuation of our approach would be the generalization of the results to 
the models perturbed far away from the simple critical curve ( |3.5| ). As we showed in section 2, 
the instanton approach works well for any two matrix model, for a general potential. Thus our 
methods are directly applicable to study the non-perturbative effects in the flows between different 
minimal models. 

In conclusion, the non-perturbative effects in the large N matrix models seem to find their 
most natural and universal interpretation in the eigenvalue instanton method studied here. Its 
application range should be much wider then the methods based on orthogonal polynomials and 
could be useful everywhere where the model can be described in the quasi-classical limit by its 
algebraic curve. 
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